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ANALYTICAL LIFETIME STUDIES OF A CLOSE-LUNAR SATELLITE 


By William R. Wells 
Langley Research Center 


SUMMARY 


An analytical study has been made to determine the influence of the third 
zonal harmonic of the moon on the lifetime of close-lunar satellites. The 
orbits considered are applicable to a close-lunar photographic mission and for 
a mission to determine the magnitude of the coefficient of the third zonal 
harmonic of the moon Jjo* Th e lunar model used in the analysis has oblate- 
ness, equatorial ellipticity, and a nonsymmetric distribution of mass in the 
northern and southern hemispheresj that is, the gravitational potential func- 
tion includes terms involving the coefficients of the zonal and sectorial, 
harmonics J 20 > ^ 22 > an & J30 in addition to the point-mass contribution. 

Target-site selection for a photographic mission is considered for orbits 
with initial arguments of pericentron corresponding to optimum lifetimes based 
on a given magnitude of Jjq (independent of its sign). With this constraint 
imposed on the initial argument of pericentron, the effect of inclination on 
the target sites and rate of change of pericentron altitude is considered. 

This phase of the analysis is considered by using several values for the mag- 
nitude of J50 (Pericentron is defined as the point in the orbit of one body 

about another at which the least separation of the two bodies occurs.) 


INTRODUCTION 


In order to design proper orbital elements for lunar missions having 
close-approach orbits, it is necessary to obtain an estimation of possible 
maximum perturbations on the satellite. The major perturbations arise from 
the disturbing influence of the earth and sun and from the higher order har- 
monics of the lunar gravitational potential function. A principal perturba- 
tion on a close-lunar satellite is due to the mass anomalies which are 
expressed mathematically in the gravitational potential function in the form 
of the coefficients of the higher order harmonics such as J20 ^ 22 ’ ant * J30 

The harmonic with coefficient Jjq causes a first-order change in the eccen- 
tricity of the orbit and consequently in its pericentron altitude. The har- 
monics with coefficients J20 and J22 cause no change and second-order 
changes, respectively, in the eccentricity, and for this reason only the effect 
of the third zonal harmonic on the satellite lifetime will be considered. 



Until a lunar satellite is established there will be no data available 
for making analytical calculations of the J^o coefficient for the moon. The 

situation is analogous to that for the earth, for which the pear-shape effect 
was discovered after analysis of tracking data from the first earth satellites. 
The best that can be done at present is to make estimates of the J^o value 

for the moon based on knowledge of corresponding values for the earth. These 
approximations may turn out to be rather gross but they do provide a starting 
point . 

It has been suggested in reference 1 that estimates of the coefficients 
of the higher harmonics of the lunar gravitational potential can be based on 
the assumption that the moon and earth can support equal stresses. For a 
smaller and less massive body, a given stress implies greater gravitational 
anomalies than for a larger body. A simple dimensional analysis on the stress 
and the equal stress assumption indicate that corresponding coefficients of 
the moon can be some 36 times those for the earth. The maximum value of J-jq 

for the earth has been estimated from studies of earth- satellite tracking data 
to be about - 2.6 x 10“6. This means that the magnitude of J^q for the moon 
may be as large as 9-3 X 10"^ if the equal stress assumption is true. 

In the absence of a Jjq term, close satellites of the moon would have 

essentially an indefinite lifetime (ignoring solar and earth perturbations). 

The presence of a J^q term causes long periodic perturbations (period of 

about a year) which will limit the lifetime of these satellites; the amount of 
limitation depends on the initial conditions under which the orbit is estab- 
lished. Since long lifetimes are generally desirable for most missions, it is 
natural to look for situations for which long lifetimes are obtained even 
though the magnitude and sign of are not known. At the same time, since 

an analysis to determine the lunar gravitational field is part of any mission, 
the orbits should be designed so that the value of J^ 0 can he estimated from 

tracking data. It will be shown in this analysis that these two features can 
be obtained simultaneously by a proper choice of the initial value of the 
argument of pericentron which is the essential parameter in the analysis. 


SYMBOLS 


a semimajor axis of Keplerian orbit, km 

* 

E time rate of change of a general element 

e eccentricity of a Keplerian orbit 


^30^m sin 2 i - e 2 cos^i 
2J 20 a(l - e 2 ) e sin i 
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G 


universal gravitational constant 


i orbit inclination referred to lunar equator, deg 

J 20 coefficient of second zonal harmonic of lunar gravitational potential 

function, taken as 2.073 X 10“^ 

J 22 coefficient of second sectorial harmonic of lunar gravitational 

potential function, taken as -2.03 X 10~5 

J-^O coefficient of third zonal harmonic of lunar gravitational potential 

function 

Jnm > K-rrm general coefficients in lunar gravitational potential function 
M mean anomaly, deg 

n mean motion of lunar satellite in Keplerian orbit, rad/day 

nj_ mean motion of moon about its polar axis, 0.23 rad/day 

Pnm associated Legendre function, Pnm( s:i - n 0) 

r* position vector from center of moon to satellite, km 

r distance from center of moon to satellite, km 

rp distance from center of moon to satellite at pericentron, km 

Arp increment of change in pericentron distance, km 

R disturbing function due to mass anomalies, km^/day^ 

Rm mean radius of moon, 1738.1 km 

t time, days 

T satellite lifetime, days 

U lunar gravitational potential function, km^/day^ 

v true anomaly, deg 

x,y, z Cartesian coordinates of satellite referred to moon-centered axis 
system, directions fixed in space, km 

x^y^z 1 Cartesian coordinates of an element of mass, referred to moon- 
centered axis system, km 

X^ m Hansen coefficient 

7 vernal equinox 

A operation denoting incremental change 
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Kronecker delta function 


0 longitude of satellite measured in equatorial plane of moon from the 

mean earth-moon line, positive eastward, deg 

ji product of gravitational constant and mass of moon, 3*6601 x lCr - ^, 

km^/ day^ 


a longitude of sun relative to mean earth-moon line, deg 

0 latitude of satellite measured from the lunar equator, positive 

northward, deg 

co argument of pericentron, deg 


“s = 


a2(l - e2) ; 


-(l - j- sin^i), secular rate of argument of pericentron, rad/day 


(^^cOm 


initial argument of pericentron corresponding to optimum lifetime, deg 


Q, longitude, from vernal equinox, of ascending node of Keplerian orbit, 

deg 

ft 1 longitude of ascending node measured from the mean earth-moon line 

positive eastward, deg 


Subscripts: 

E earth 

0 initial value 

1 short-period terms 

2 long-period terms (half the rotational period of the moon about its 

polar axis) 

5 long-period terms (period of the line of apsides of the lunar orbit) 

k secular terms 

A dot over a symbol represents differentiation with respect to time. 

A bar over a quantity in parentheses denotes the mean value of that quantity. 


GENERAL CONSIDERATIONS 


Lunar Model and First-Order Perturbation Analysis 

The mutual attraction between a body with finite size and a body con- 
sidered as a point mass can be studied by analyzing the motion of the point 


mass relative to the finite body. This motion is expressed mathematically as 


^(x,y,z) = (l) 

where r? is the position vector of the satellite or point mass relative to the 
mass center of the moon or finite mass and is the mutual force of attrac- 
tion, per unit of satellite mass, between these two bodies. 

The force field I? is conservative since it is due to gravitational 
attraction only and, as such, can be represented as the gradient of a scalar 
U called the gravitational potential function, that is. 


F = W (2) 

or 

r = W (3) 

The gravitational potential function per unit satellite mass for each ele- 
ment of the attracting body is directly proportional to the element of mass and 
the universal gravitational constant and inversely proportional to the distance 
separating this element of mass and the satellite. The potential function per 
unit of satellite mass for the totality of elements of the finite body is then 
expressed as 


u * ° JF F — 2 — w* (k) 

m ftx - x ') 2 + (y - y *) 2 + U - z') 2 J 

where M in equation (4) is the mass of the attracting body. 

The potential function can be expressed in terms of the latitude and lon- 
gitude of the satellite relative to coordinates on the surface of the moon and 
the distance of the satellite from the center of mass of the moon. This expres- 
sion can be made by recognizing that the gravitational potential satisfies 
Laplace's equation for points exterior to the moon's surface and by solving 
this equation for U in terms of spherical harmonics (ref. 2). Laplace's 
equation is 

= 0 (5) 

The solution of equation ( 5 ) in terms of the notation of the present paper is 


U(r,9,0) 



00 n ^ 

1 m « 

n =1 ttj=0 


cos m9 + K nm sin m0)P nm (sin 0) 


( 6 ) 
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The solution of Laplace's equation as given by equation (6) is given in refer- 
ences 2 and 3* In the present analysis it is assumed that the potential func- 
tion given by equation (6) can be approximated by the terms having coeffi- 
cients J 20 > ^ 22 > ant ^ ^30* These harmonics account for whatever oblateness, 

equatorial ellipticity, and pear-shape is present. This assumption gives the 
following expression for the potential function: 


u(r,e,0) = H 


1 + J 2 o( 1 - 3 sin 2 )*) 


2r 


- J22 c °s 2 0 cos 20 J 2Q sin 0^5 sin 2 0 - 3) 

r 2 2r3 s 


(7) 


The first term on the right-hand side of equation (7) is the potential due 
to a point mass; the second term accounts for the oblateness of the moon; the 
third term accounts for the ellipticity of the equator, and the fourth term 
accounts for a nonsymmetric distribution of mass in the northern and southern 
hemispheres (pear-shape effect). 


The first-order perturbations of the elements (excluding mean anomaly M) 
of the close-lunar orbit on the values obtained for a point mass are developed 
in appendix A and are given here for reference and include both long-period 
effects and the secular effects: 


Aa = 0 


Ae = - — sin i(sin co - sin o> 0 ) 


2aJ 


20 


Ad = 


3nJ2Q^m 

a 2 (l - e 2 ) 


2! - 1 sin% ) t + 


3nJ 22 Rm 


2nja 2 (l - e 2 ) 


1 - — sin 2 i)(sin 2 ft' 


L )(si ; 


(8a) 

(8b) 


sin 2 ft, 


;) 


J 30 R m 


2aJ 20 (l - e 2 ) 


sin^i - e^cos^i 
e sin i 


(cos CO - cos C0q) 


(8c) 


e cos i 


3 nJ sin i f x 

Ai = -(cos 2 ^ T - cos 257 ^) + ~(sin co - sin co Q ) ( 8 d) 

2 a^n^_ (l - e 2 ) 2 J_ 20 a ( 1 " e / 
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(8e) 


) 


J^O R m e cos i 

J (COS O) - COS C0 o ) 


The elements a, e, and i that appear in the coefficients above, strictly 
speaking, are mean values as discussed in appendix A. 

Most of the angular quantities in the preceding equations are illustrated 
in figure 1. In the following discussion perturbations of differing durations 
are discussed. Perturbing terms involving true anomaly v are said to cause 
short-period effects since these effects are periodic with a period equal to 
that of the satellite in its orbit. Perturbing terms containing the longitude 
of nodes relative to the mean earth-moon line ft 1 are said to cause long- 
period effects which have a period approximately equal to half the period of the 
moon’s rotation about its polar axis. Terms involving the argument of peri- 
centron 0) are said to cause long-period effects which have a period equal to 
the period of the line of apsides of the lunar satellite orbit. Perturbing 
terms which are independent of the angles v, Q x , and co are referred to as 
secular effects. These effects are permanent and, as such, have an infinite 
period. In the expressions for the variations in the elements given by equa- 
tions (8), the short-period effects have been averaged out leaving only the 
secular and both long-period effects. The changes in the elements due to short- 
period effects are, in general, much smaller than the long-period and secular 
effects and will be considered in a later section of this paper. 

As indicated in equations (8), the mass anomalies cause no first-order 
perturbations in the semimajor axis. The arguments of pericentron and longitude 
of node, however, have perturbations due to all three effects; that is, 
oblateness, equatorial ellipticity, and pear-shape. The effect of the oblate- 
ness (jgO term) on these two elements is a secular change. The effects of 

equatorial ellipticity and pear-shape ^22 an< ^ J 30 terms ) take the form of 

long-period changes. The eccentricity 
has only a long-period change caused by 
the pear- shape effect. The change in 
inclination is due to both long-period 
effects - equatorial ellipticity and 
pear- shape. 

Lunar equatorial plane 

Variation in Pericentron Altitude 

r 

Memal equinox 

In a study of close-lunar satel- 
lite lifetimes, the main concern is the 
manner in which the pericentron altitude 

changes with time. For an orbit in Figure 1.- Illustration of pertinent angles. 




JnJgQ^n cos i . 5^22^m cos ^ 

2a 2 (l - e 2 ) 2 2n^a 2 (l - e 2 )^ 
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which the s emima jor axis is constant (for instance, when mass anomalies cause 
the perturbations and with short-period effects neglected), the change in peri- 
centron altitude is given by- 


Arp = -a Ae ( 9 ) 

A first-order approximate analysis can be performed on the satellite life- 
time with the change in eccentricity, given by equation (8b), substituted into 
equation ( 9 ) which gives 


Arp = i R m sin i(sin to - sin co 0 ) 


t= - ^2. + — sin - "^ 
cb 05 


2J 20 Ar p 
iJ30 R m sin 1 


+ sin cD r 


( 10 ) 


If the expression for cko/dt obtained in appendix A Is used, the variation in 
Arp with time can be obtained. It is shown in appendix A that an approximate 
expression for cko/dt is given by 


do 

dt 


3 nJ 20 R m / 
a 2 (l - e 2 ) 2 ' 1 


l sin2i ) 


( 11 ) 


The main analysis will be made (with the exception of a comparison of lifetimes 
using first- and second-order results) by utilizing the first-order results of 
equations (10) and (ll). 


A second-order expression for Arp is developed in appendix B. This 

result is obtained by assuming cko/dt is given by the sum of the secular and 
long-period terms, that is 


cko 

dt 


3nJg 0 Rg 

a2(l - e 2) 2 



5 , 2 ^ 
t sln ^ 


1 + 




m 


sin to 


sin 2 i - 


e 2 cos 2 i 


2J 20 a 


(l - e 2 ) 


e sin i 


( 12 ) 


The result of this analysis is that Arp, including second-order terms, is 
given as 


where 


j 30 R hi s: *- n 1 .. fl + f sin co 
Ar-n = ■ log 


2fJ 


20 


ll + f sin ov 


f = 


^30^™ sin 2 i - e 2 cos 2 i 


2J20 a 




e sin i 


( 13 ) 
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Analysis of Short-Period Perturbations and Earth Effects 


The short-period variations in the elements due to the lunar oblateness 
are developed in reference k. In the notation of the present analysis, the 
short-period variations of semimajor axis and eccentricity are written as 



| sm" 


■) (t f - - 




sin% cos 2 (a) + 



v 


v 0 




i^r^sin 2 i COS + 


3 sin 2 iJ 20 
4a 2 e(l - e 2 ) 


cos 2 (<m + v) 


( 15 ) 


+ e cos(v + 2m) + ~ e cos(3v + 2m) 


(16) 


The change in the pericentron altitude due to these short-period varia- 
tions is obtained by substitution of equations ( 15 ) and (l6) into the following 
expression for ^Ar V : 


= (! - e)( Aa ) i " 


( 17 ) 


The variations in the pericentron altitude of the lunar orbit, due to the 
disturbing effect of the earth, can be computed from the following equations. 
The mean rates (over a lunar day) of argument of pericentron, eccentricity, and 
semimajor axis due to the influence of the earth are: 




= — e\]l - e^ sin^i sin 2m 

8 n v 

= 0 


(18) 

( 19 ) 

( 20 ) 
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For close-lunar satellites moderately inclined, i s Sma - L - 1 - compared 

with (jd s . If dn/dt is approximated by cb s , then equation (19) can be inte- 
grated to obtain (Ae)g as: 

(Ae) E = - e\/l - e 2 sin 2 i(cos 2x> - cos 2l u G ) (21) 

l6no> s 

Substitution of equation (21) into equation (19) gives for ^Ar^j ^ the fol- 
lowing expression: ' 

an 2 

(At e]fl - e^ sin^i(cos 2 d - cos 2x>o) (22) 

\ P / ■“ 16 nm„ 


From equation (22) it can be deduced that the maximum decrease in the peri- 
centron altitude due to earth effects is 


(Ar \J| = _ — — i. e^l - e^ sin% 

1 + 

COS &J 0 o 1 

{ pjEJ max 16 v 




(23) 


The maximum value of j^Ar^ max occurs for cu 0 = 0° or l8o°, and is 

[(^sjmax = ' ^8 * <T " ® 2 Sin2i (2i+) 


LUMP SATELLITE LIFETIME STUDY 


Lifetime Analysis 

In this analysis the lifetime of a close-lunar satellite is considered 
ended whenever the altitude of pericentron decreases more than a specified 
amount which depends on the initial pericentron altitude. For instance, the 
lifetime of an orbit with an initial altitude of pericentron of 46 kilometers 
might be considered ended whenever the decrease in rp is more than 36 kilo- 
meters, that is, when Arp = -3 6 kilometers. This criterion on lifetime would 

allow a 10-kilometer variation in pericentron altitude which may be present as 
a result of a combination of effects neglected in the analysis such as second- 
order effects in u>, short-period perturbations, and perturbations due to the 
earth. 

The first-order analysis on the lifetime of the satellite is obtained from 
a solution of equations (10) and (ll) which are repeated here. 
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Ar = 

P 2 J 20 


R m sin if sin o> - 


i ( £ 


sin u> 0 ) 


du> 

dt 


***& -(! - \ 

a 2 (l - e 2 ) 2 ' 4 


In all examples the orbit used (unless otherwise noted) is a direct orbit for 
which the pericentron and apocentron altitudes are 46 and 925 kilometers, 
respectively. (Apocentron is defined as the point in the orbit of one body 
about another at which the greatest separation of the two bodies occurs.) The 
inclination is 21 ° and the initial location of the node is 50° west of the mean 
earth-moon line = - 50 °) • This situation corresponds to a typical earth- 

moon transfer trajectory for which the transit time is about 80 hours. The 
magnitude of J^q is assumed to be 9.3 x 10 “ 5 . a parametric study on both 
inclination and J-jq is performed. 


The variation in rp with time as given by equations ( 10 ) and (ll) is 
illustrated in figure 2 for a negative sign on J30 T * le efTect of some ini- 
tial values of the argument of pericentron is illustrated. A plus sign for 
J30 simply causes a change of sign on the values presented in figure 2 . These 
curves indicate that the pericentron altitude can either increase or decrease 
initially depending on the values of cjq 0 and J30 that are used. 

A solution for the satellite lifetime T from equations ( 10 ) and (ll) is 


illustrated in figure 3 for a range of 

esting features are present in this 
figure. With J ?0 = - 9-3 X 10 - 5 , a 
minimum lifetime of about 2 weeks 
occurs for cjd q near 0° and infinite 

lifetimes occur for the range 
4 - 8 ° < <x> 0 < 132 °. A similar situa- 
tion exists for J30 = 9*3 X 10 “^ : a 

minimum lifetime of 2 weeks occurs 
for <x> 0 near l8o°, and an infinite 
lifetime occurs for 228° < <x> 0 < 312°. 

The explanation for the shift 
from a finite to an infinite lifetime 
can be obtained from figure 4 -. For 
o) 0 < 48 °, for instance cjd q = 30°, the 
curve intersects the line marked 
Arp = -36 kilometers after about 
16 days; hence, a l6-day lifetime is 
shown in figure 3. The curve for 
o^o = 48 ° is tangent to the line 
marked Arp = -36 kilometers after 


03 o from 0 ° to 360°. Several inter- 



Figure Z- Variation of Ar p with time for several initial values of the 
argument of pericentron. 
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about 35 days so that this 
time serves as a limiting 
value of the finite lifetimes 
for cd 0 < 48°. For 
48° < ojq < 132°, for instance 
0) o = 60°, the curve never 
intersects the line marked 
Arp = -36 kilometers; hence, 

the infinite lifetime. 


The estimate of the mag- 
nitude of Jjq of 9*3 X 10“ ^ 

provides a possible basis for 
determination of the magni- 
tude of the pericentron alti- 
tude variations. However, 
the sign of J^q is very 

critical to the actual life- 
time of the satellite. Since 
the sign is unknown, it is 
natural to inquire into the 
possible existence of initial 
conditions on the satellite 
that result in "optimum" or favorable lifetimes for a given magnitude of J 

regardless of its sign. Reference to figure 3 indicates that for 
| J30 1 = 9*3 X 10“5 there is an initial value of the argument of pericentron 

of approximately 48° for which 
the lifetime is infinite if 
JjO < 0 an d is favorable 

(T = 95 days) if Jjq > 0. The 

finite value is called the "opti- 
mum lifetime." A satellite in an 
orbit having this initial value 
of argument of pericentron will 
have at least a lifetime of 
95 days . A similar situation 
exists at co 0 = 228°. For the 
remainder of the analysis, when 
optimum lifetimes are considered, 
the initial value of argument of 
pericentron less than 90° will be 
used and referred to as (co 0 )m* 

It should be noted that 
(<n 0 )m = ^8° was for the special 

Figure 4.- Illustration of infinite and finite lifetimes concept 



Time, days 


280r 


46-925 km orbit 

i = ? 1° 



Figure 3.- Satellite lifetime as a function of initial values of argument of 
pericentron. 
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is dependent on 


case of i = 21° and |J-^q| ” 9*3 X ICf^. In general (to 0 ) m 
both inclination and I J^q [• 


Photographic Mission 

A lunar photographic mission to he performed by a satellite in a close- 
lunar orbit should be designed so that the orbit has a sufficiently low peri- 
centron altitude to provide good photographic resolution. At the same time, 
the pericentron altitude must be sufficiently high so that the satellite life- 
time is not ended before the photographic mission is completed. Other charac- 
teristics of the orbit such as eccentricity, inclination, and initial argument 
of pericentron must also be chosen properly since they affect the lifetime of 
the satellite as veil as the photographic coverage of the lunar surface. The 
value of Initial argument of pericentron greatly affects the lifetime of the 
satellite and the area available for the photographic mission. For this reason 
it is chosen to correspond to the value that results in the optimum lifetime, 
discussed previously. 

Another factor to be considered in a photographic mission is that of 
favorable lighting conditions. The lighting conditions on a region are con- 
sidered favorable for photography whenever the sun rays are at about 60 ° to 
the local vertical. This situation is illustrated in figure 5 . As shown in 
figure 5 , the Initial argument of pericentron should be situated near the point 
where the trace of the lunar satellite plane intersects the 60 ° circle centered 
about a line in the direction of the sun. These conditions will insure favor- 
able lighting for the photographic phase of the mission. 


N 



Figured- Illustration of initial angular relationship for photographic 
mission. 
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Effect of Orbit Inclination on Target Site Selections 
and Lifetimes of Close-Lunar Satellites 

Utilization of the optimum lifetime (previously defined) in a lunar mis- 
sion automatically constrains any operation dependent on the initial value of 
the argument of pericentron. In the case of a close-lunar photographic mission 
in which the pictures are to be taken at pericentron, this constraint takes the 
form of a restriction in surface areas available for the experiment. This con- 
straint, however, can be made less severe by allowing some freedom in the selec- 
tion of the inclination for the mission. The effect of inclination on changes 
in pericentron altitude, optimum lifetimes, initial latitudes and longitudes of 
pericentron location, and subsequent changes with time, initial solar positions, 
and initial values of (<u 0 ) m is considered in this section. In all cases, the 

orbits considered are direct and have pericentron and apocentron altitudes of 
46 and 925 kilometers and = -50°. The values of (o> 0 ) m were determined 

from the assumption that the life of the satellite ended whenever the peri- 
centron altitude decreased more than 36 kilometers. 




80- 



i, deg 


Figure 6.- Initial values of argument of pericentron. 


The condition needed to obtain an 
expression for (oo 0 )m can be deduced from 

the definition of and from refer- 

ence to figure 4. It can be seen from this 
figure that o>o = (o) 0 ) m 'whenever the condi- 
tions ^^Arpj = 0 and Arp = -36 kilometers 

occur together. If this operation is per- 
formed on equation (lO) and if the value 
of ATp = -36 kilometers is used in equa- 
tion (lO) it is found that 



The dependence of (^o)m on ^ as 

determined from equation ( 25 ) is shown in 
figure 6. The maximum value of (o) 0 ) m for 

J^O = -9*3 X 10“ 5 is about 65 ° corre- 

ponding to i = 90°. The minimum value is 
0° and occurs for i = 5°. For inclina- 
tions less than about 5°, the lifetime of 
the satellite is unaffected by the presence 
of a J^o term of this magnitude. Fig- 
ure 6 also illustrates the variation of ini- 
tial latitude and longitude of pericentron 
location with inclination. Note that these 
initial longitudes and latitudes are 
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essentially restricted to the northwest (cartographic) quadrant of the moon. 
Other areas are possible if a different direction in orbit is considered. 

The change in pericentron altitude with time, for various inclinations, 
is shown in figure 7 • These curves were generated with gd 0 = (^o)m an( ^- 

J 30 = -9«3 X 10“5. Note that Arp is never less than -3 6 kilometers since, 
according to figure 3 , these values of J^q and od q = cause an infi- 

nite lifetime. However, if a value of J^q = 9-3 X 10“5 were used, the curve 
in figure 7 "would be inverted and, for o ) 0 = ( CD o)m^ an infinite lifetime would 

not occur. (See fig. 3 for the case of i = 21°. ) Instead, the optimum life- 
times which occur are, for example, 104 days for i = 15 °, and 124 days for 
i = 10 °. 


Figure 8 presents the variation in the optimum lifetime with inclination. 
Inclinations of less than about 5° result in an infinite lifetime. Inclina- 
tions greater than 5° result in finite optimum lifetimes. The minimum of the 
optimum lifetime, about 95 days, occurs for an inclination near 25 °. 


The initial rate of change of the peri- 
centron altitude for go 0 = (oo 0 ) m and 

J30 = - 9-3 X 10~5 i s also shown in figure 8 , 

For 5° < i < 60 °, the pericentron decreases 
initially. For 5° < i < 15°.> this decrease 
continues until about 45 days have passed, as 
can be seen from figure 7 - The initial rate 
of change in pericentron altitude for 
i = 10 ° and 15 ° is - 1*3 an d 
-1*7 kilometers/day, respectively. Figure 7 
indicates that these rates remain approxi- 
mately constant for the first 10 days in 
orbit. These rates indicate that in addition 
to the provision of optimum lifetimes, the 




Figure 7.- Variation of pericentron altitude with time. 


Figure 8.- Initial values of rp and 0 and optimum 
lifetimes. 
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condition aa 0 = (^o)m causes changes in the pericentron altitude that could 

reveal, through an analysis of the tracking data, an early determination of the 
sign of J^o an< ^ a first estimate of its magnitude. 

The initial solar position cr 0 is also shown in figure 8. An initial 
longitude of node Oq = ~ 50 ° and solar incidence angle of 60 ° to the local 
vertical were used in the determination of the solar positions. The value of 
a Q varies from 10° east to about 57° east for inclinations between 5° and 30 °. 
These angles correspond to about 1 and days before full moon. Reference to 

figure 5 indicates that it is also possible to obtain favorable lighting con- 
ditions for these initial elements after full moon. The values of cr 0 in fig- 
ure 8 and the values of 0 O and 0 O in figure 6 were computed by use of equa- 
tions (C2) and (C3) from appendix C. 


The pericentron locations, latitude and longitude, are shown in figure 9 
as a function of time after orbit establishment. Figure 9 shows that both 0 
and 0 essentially change linearly with time; the rate of change in 0 is 



Figure 9.- Variation of pericentron location with time 
after orbit establishment. 


much less than that in 0. A corre- 
sponding plot of surface coverage is 
shown in figure 10. Note that the con- 
straint of optimum lifetime has made 
unavailable, to a photographic mission, 
the greater part of regions in the 
immediate vicinity of the lunar equator. 
A means of including this region , how- 
ever, could be obtained by such opera- 
tions as a delay in orbit or a simple 
plane-change technique. 

Effect of the Magnitude of J^q on 

the Lifetime of a Close-Lunar 
Satellite 

The previous analysis was performed 
by using the value Jjq 5 -9*3 X 10~5. 
Since this value is only a first esti- 
mate of J-^q, it is of interest to 
investigate some of the previous results 
when other values of J^q are assumed. 

The effect of J^q on such parameters 
as (^ p ) max j ( a) o)m^ optimum lifetimes, 
(*p)o, an< i pericentron location are 

shown in figures 11 to 13- The values 
of Jjq used are the suggested value 
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L- 64-1018 

Figure 10.- Trace on Moon of argument of pericentron with time. 
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^ r P^max, km 


of 9.3 X 10 - 5, twice this value, one-half this value, and one-fourth this 
value. In all cases cd 0 = (<^>o)m* 

Figure 11 indicates that there are certain values of inclination, 
depending on the value of J^ 0 , for which (u> 0 ) m is zeroj values of i below 

these correspond to indefinite lifetimes for the satellite. In general, the 
inclination for (<x> 0 ) m = 0 decreases with an increase in magnitude of J jq . 

Equatorial orbits (i = 0°) were not considered in the analysis since one 
requirement of the mission is the determination of a value for J30* Figure 13 

gives the location of the pericentron for the first week after establishment 
of the lunar satellite corresponding to three values of J 30 and i = 15°. 

For a variation in Jjq from -4.65 x 10" 5 to -1 8.6 x 10"5, a change in (CD 0 ) m 

occurs which causes the initial latitude of pericentron to change from about 
4° to about 11°. The corresponding changes in initial longitude of the peri- 
centron are from -34° to 4°. 






Figure 11.- Variation of (Ar p ) max and (w 0 ) m with inclination. Figure IZ- Variation of (f p ) 0 and optimum lifetime 

with inclination. 
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INFLUENCE OF SECOND- ORDER EFFECTS, SHORT-PERIOD PERTURBATIONS, 

AND EARTH EFFECTS 


Second- Order Effects 

The second-order effects of J^q on the lifetime of the lunar satellite 

are determined from equations (12) to (l4). These effects are dependent on the 
magnitude of the correction function f given by equation (l4). This function 
is seen to be dependent on a, e, i, Jjqj an & R m* For a given orbit 

size and assumed lunar gravitational potential, the value of f varies with i 
as shown in figure 14 . In the vicinity of i = 10°, second-order effects are 
negligible but may become important at lower and higher inclinations. The 
effect of this correction on the lifetime of a satellite is shown in figure 15 . 
In general, the lifetime decreases with the greatest effect (a 20-day differ- 
ence) near co 0 = l80°. The band of values of cd 0 for which the satellite has 
an infinite lifetime also narrows. 




Time, days 


Short-Period Perturbations 
and Earth Effects 

Some results for short-period 
change in pericentron altitude are 
computed from equations ( 15 ) to ( 17 ) 
and are shown plotted in figure 1 6 . 

The change in pericentron altitude 
from the time the satellite is at some 



Figure 13.- Variation of pericentron location with time for var- 
ious values of J3Q. 


Figure 14.- Maximum effect of correction to u) 
with inclination. 
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initial point in the orbit v = v Q to the time it is at pericentron v = 0° 
is shown for all possible initial true anomalies. Two initial arguments of 
pericentron are considered, a.) 0 = 0° and od 0 = 40° for each of two orbits: 

one having a pericentron of 46 kilometers and an apocentron of 925 kilometers, 
the other having a pericentron of 46 kilometers and an apocentron of 1850 kilo- 
meters. The largest decrease in pericentron altitude (less than 0.6 kilometer 
for these orbits) occurs for an initial location of the satellite in the vicin- 



deg 


ity of apocentron 
(v Q = l80°). These vari- 
ations in pericentron 
altitude are small com- 
pared with the long- 
period variations obtained 
by assuming the suggested 
magnitude of 9*3 X 10 “ 5 
for J^q. 

The effect of the 
earth on the lifetime of 
a close-lunar orbit can be 
determined from equa- 
tions (l 8 ) to (24). In 
figure 17 , the variation 
in the pericentron alti- 
tude with time due to the 
earth effects is shown. 


Figure 15.- Comparison of lifetimes for first- and second-order analysis. 



v 0 , deg 


Figure 16.- Short periodic variation in pericentron altitude with time. 


The maximum decrease in 
the pericentron altitude 
for these initial argu- 
ments of pericentron is 
about 5*0 kilometers. The 
maximum value of 


[(%)£] max is given by 

equation (24) as a func- 
tion of inclination for a 


given orbit size. This 
value is shown plotted 
in figure l 8 against i 
for a 46-925- and a 
46 - 1850 -kilometer orbit. 
For these orbits it can 


be seen that the maximum 


value of [(Arp^nax 


lies in the range from 0 
to 70 kilometers for 
inclination between 0° and 
50°. For the 21° inclined 
orbit the average values 
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Figure 17.- Variation of pericentron altitude with time due to earth 
effects. 
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these orbits are about 6 and 14 kilo- 
meters. These results indicate that Figure 18. - Maximum change in pericentron altitude due to earth 

both the short-period and earth per- 
turbations on the pericentron altitude 
of low inclined close-lunar satellite 

orbits are small compared with the previously discussed long-period effects but 
can become significant for more inclined and larger orbits. 


CONCLUDING REMARKS 


An analytical lifetime study of close-lunar satellites has been performed 
for a mathematical model of the moon which has oblateness, equatorial ellip- 
ticity, and a nonsymmetric distribution of mass in the northern and southern 
hemispheres. It is shown that a value of the coefficient of the third zonal 
harmonic J^q for the moon of approximately times that for the earth can 
cause large changes in the pericentron altitude resulting, in many instances, 
in short lifetimes of close-lunar satellites. On the other hand, it is shown 
that there are conditions for which the lifetime is indefinite for a favorable 
sign of Jjo- 

In general, an arbitrary initial argument of pericentron either results in 
a long lifetime if Jjq has one sign or a short lifetime if is of the 

opposite sign. However, it is shown that, associated with each orbit, there 
are two initial arguments of pericentron, (^o)m and (u> 0 )m + l 80 °, for which 

the lifetime is relatively long independent of the sign of J-^q. Utilization 

of these f, optimum lifetimes" in a lunar mission automatically constrains any 
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operation dependent on the initial value of the argument of pericentron. It is 
found that, in the case of a close-lunar photographic mission in which pictures 
are to be taken at pericentron, this constraint takes the form of a restriction 
on the lunar-surface coverage available for the experiment . However, this 
restriction can be lifted somewhat by allowing some freedom in the choice of 
orbit inclination for the mission. Under this restriction, the region at the 
center of the lunar disk is unavailable for immediate photographic coverage in 
a typical close-lunar orbit. 

It is shown that the constraint of optimum lifetimes allows a great deal 
of freedom in selecting photographic surface areas and a means of obtaining 
data on the lunar gravitational field from a tracking analysis. The magnitude 
of the time rate of change of the pericentron altitude should be sufficiently 
large for this condition to obtain a first estimate of the magnitude of Jjq 
as well as its sign. 

The general results of the analysis are based on satellite lifetimes as 
determined from a secular change in the argument of pericentron. The lifetimes 
as determined by a more accurate representation of the change in argument of 
pericentron, a secular plus a long-period change, is considered briefly. The 
result of the calculations indicated that the more accurate analysis results in 
shorter lifetimes. 

Results of an analysis of short-period perturbations and earth effects 
indicated that the lunar satellite lifetime is not affected, to a great extent, 
by these factors whenever the orbits are moderately inclined and are not too 
large . 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 15 , 1965. 
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APPENDIX A 


PERTURBATION OF THE ELEMENTS 


Analyses similar to the following ones are performed in numerous papers. 
(See, for example, ref. k.) The intent here, however, is to give some of the 
details not present in the other developments. 


The gravitational potential function used to represent the moon in this 
analysis is 


u(r,e, 0 ) 


1 + J2o( 1 - 3 sin 2 0) 


2r 


Jop cos 2 0 cos 20 - J^o sin 0(5 sin^0 - 3 ) 
r d 2x-> 


(Al) 


The disturbing function due to the fact that the moon has mass anomalies is 


= u - i £ = 


2 r; 


^|£(l - 3 Sin 2 0) 


- 5J22 cos 2 0 cos 29 - ^ J 30 s i n 0(5 sin 2 0 - 3 ) 


(A2) 


The disturbing function R can be expressed in terms of the elements of the 
orbit from a consideration of the following sketch; 
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Application of spherical trigonometry to the geometry presented in the sketch 
gives 


cos 9 = — t[cos(cd + v)cos ft' - cos i sin ft'sin(cD + v)l 
cos <p'— - 1 


(A3) 


From this expression the term cos 2 0 cos 29 in equation (A2) can be expressed 


as 


cos 2 0 cos 29 = — (l + cos 2 i)cos 2T2*cos 2(cd + v) 


2 _ 

+ sin 2 i cos 20,' - 2 cos i sin 2ft' sin 2(co + v) 


) 


(A4) 


In terms of the elements of the orbit, R can be written as 
R = - 3 sin 2 i + 3 sin 2 i cos 2 (cd + vj] 

- ^ ^ |(l + cos 2 i)cos 2^'cos 2(cd + v) + sin 2 i cos 20' 

- 2 cos i sin 2ft' sin 2(o> + v) 


- 5(5 sin 2 i - U)sin(cu + v) 


- 5 sin 2 i sin 3 (to + v 


)j sin ij 


(A5) 


To eliminate the short periodic perturbations (on the order of one orbital 
period), integrate R with respect to M from 0 to 2 n assuming that a, e, 
i, cd, and O' are constant. This procedure gives, the first order in J20> 
J 22 > a nd J 30 the following: 

S = r 2 ” ^(1 - e2)- l/2 a, 

J q 2tC J q \ a / 

, - e^ /2 {x - | sin2i) - ^(l - e2)- 3/2 sin2i 


COS 20' 


2a 30 


(l - e 2 )" 5//2 sin i(: 

- ^ sin 2 i)sin co 


4 / J 


(A 6) 


The integration indicated in equation (A6) can be determined immediately if use 
is made of the following results obtained from reference 5 • 
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-yU y HI 

A o 


cos mv dM 


- 1 f 0 ”(ir 




sin mv dv = 0 


(A7) 




where p is a summation index. 


To find the variation in the elements due to this averaged disturbing 
function TZ, the following LaGrange equations are used (ref. 6): 


da = _2_ 5R 
dt na dM 




de 

dt 


ckjp 

dt 


1 - e" 


£ Sr _ )Jl - e g dR 


na^e 


2 e 


na^e 


cos 1 dR \|l - e 2 dR 




na^Vl 


sin i ^ 


na e 


de 


cos i 


dR 


dR 


di _ _ 

^ na 2 \/l - e 2 sin i dco na 2 sin l\jl - e 2 d£2 


dii 

dt 


na‘ 




L dR 

e 2 sin i 


J 


Substitution of equation (a 6) into equations (A8) gives 


(AS) 


2* = 0 
dt 


(A9a) 


j£.-^ J 3° Ktnl ( 1 -2. ln %) 

<it 2a^ A o\ 2 ^ ^ ' 


COS O) 


(l - e 2 )‘ 


(A9b) 
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diD 3nJ 20 R S A _ 5 sin 2^ + 5nJ 22 R m 


dt a2(l. e 2) 2V 4 


S (l - * 2 )‘ 


1 - sin 2 i]cos 22 ' 


2 \ 2 


^ nJ ^° Rm - Sln ^ |ie 2 - (35e 2 + i+)sin 2 i + ^(je 2 + l) sirAi] 

8a^(l - e 2 )\ sin i 


(A9c) 


di 

dt 


3nJ pp H^ sin i 3^30 R m e cos i 

a 2 (l ' - e 2 ) J 


sin 22 + 


2s? (l - e 2 f 


^1 - ^ Bin 2 lIos a) ( A9 ( i) 


&2 

dt 


P P 

3nJ 2 QR^ cos i 3 n J22 R m cos 1 

2a 2 (l - a 2 ) 2 a 2 (l - e 2 f 


cos 22 ' 


+ — 3nJ ^£ — cot l(l - ~ sin 2 i)sin m (A9e) 

2a5(l-e 2 ) 3 V 4 ' 

From equation (A9) it is seen that du)/dt and dQ/dt have secular and long- 
period changes. All the rates of change of the elements except da/dt and 
de/dt have long -period changes (about 2 weeks) and all the rates of change of 
the elements except da/dt have a longer periodic change (period of the satel- 
lite orbit line of apsides). Some idea of the nature of the long -period change 
(long- period changes due to gd) in eccentricity and inclination can be obtained 
from equations (A9b) and (A9d) . (in the long -period analysis it will be 
assumed that the effects which are periodic in Q* have been averaged out.) If 
for the moment all quantities appearing in equations (A9b) and (A9d) are con- 
sidered constant except the arguments of pericentron and if doo/dt is approxi- 
mated by its secular value, it is easy to integrate these expressions to obtain 
long-period changes in eccentricity and inclination. To account for this 
second long-period variation in eccentricity and inclination in the integration 
of equations (A9), the rates given by equations (A9) will be expanded in a 
Taylor f s series about mean values (that is, over the period of rotation of the 
line of apsides) of inclination and eccentricity. If it is assumed that the 
values of e and i are not greatly different from the mean values, then only 

two terms of the expansion need be considered, that is, if — = E(a, e, i,Go) 

.dt 

represents any of the rates given in equations (A9), then E can be expanded 
about e and 1 to get the approximate or first-order relationship: 


>• 

E(a, e, i,u>) = E(a, e, i,co) + — 

oe 


_ _ Be + M 

e.i 8i 


_ , Si 


e, 1 


(A10) 
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where 


and 


8i = i - 
6e = e - 




(All) 


(A12) 


^Note, a = a since ^ = 0.^ Approximate expressions for 8e and Si can 

be obtained by integration of de/dt and di/dt as given by equations (A9b) 
and (A9d) with the assumption that cto/dt is approximately equal to only its 
secular value, that is, by the first term in equation (A9c) and that the values 
of e and i in these expressions are e and 1, respectively. Under these 
assumptions 


de 3^ J 30 sin i / 5 2 \ 

~ 1 - — Bin 1 COS CD 


dt 2a3 (i _ 8 2) 


-p\2\ 4 


(A13) 


and 


■3s 


di 5 nJ 30 V 5 2"\ 

* 1 - 7 - smh cos 03 


e cos i 


dt 2a? (l - X “ 


(AlC) 


From this it may be seen that 


e ~ e n - 


J 30 sln 1 


2a5(L c 


(l - e2) 


5 2-\ r 

1 - 7 - sin i ) / cos od dio 
J co 0 


-p\2 \ 4 


= - 


J 30 


o - — — sin i(sin cd - sin cd q ) 


2J 2 q 


(A15) 
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3nJ3oR^e cos i/ = o_ 

i « i 0 + — r=r — (1 - r sin21 

2e?{l - e 2 )^d) s ' ^ 


cos CD dco 


cos i 

i o + 7 rr(sin co - sin ooq) 

2J on a{l - e 2 ) 


(A16) 


Substitution of equations (AJL5) an<i (Al 6 ) into equation (A 12 ) gives for e and 
i the following approximations: 


^ J 3° R m , . . 

6 = e o + 2J^ T Sln 1 Sin 

r . J 30 Rm 1 cos 1 . 

1 = 10 - ^5 ^TTi^I sin ““ 


(AIT) 


Then 


5e = e - e = - - — — sin i sin co 
2J 20 a 


(A18) 


_ J30 cos ^ 

5 i = i - i 7 — v sin co 

2J20 a 1 - e 2 ) 


(A19) 


These values can be substituted into equation (A10) to obtain an approximate 
value of E in terms of mean values of eccentricity and inclination. The 
partial derivatives to be used in equation (AlO) are listed as follows (again, 
it is assumed that all effects which are periodic in n T have been averaged 
out ) : 


a2(l - e2) 3 


12 nJ 2 ()Rm e / 5 2 -\ 

1 - f. sin"! + 0 

of, -o\5 V k ) 


H 


(A20a) 


6 nJ 0 _R 2 e 
20 m 


" " a 2 (l - e 2 ) 5 


cos i + 0 


(J 30 ) 


( A20b ) 


be e, I 


15^20 R m 

2 a 2 (l - e 2 )‘ 


sin i cos i + 0 


(j 3 o) 


(A20c) 


8a- 
8 e e. 


(A20d) 
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be 

be 


e } I 


= 0 



M 

= °( J so) 


(A 20 e) 

(A 20 f) 


Quantities of the order of J30 in the preceding equations were not evaluated 
since, as is shown later, such terms contribute terms of the order of J 2 q^J 2 q. 

These terms are considered second order with respect to terms of the order 
of Jjo» 

Substitution of the results from equations (A 20 ) into equation (A 10 ) 
gives the following secular and long-period rates: 


cks 

dt 




i 2 (l - e 2) 2 


1-2- sin 2 ! 
4 


. 2t -2 2t / = \ 

sin x - e cos xf ± _ 5 sin 2 j) sln ^ + 0 

k I 


an 

dt 


2a3(l - e 2 ) 3 e sln 1 

3 nJ 2 Q R^ cos i ^ 3 nJ^ 0 R^e cos l(l - £ sin 2 I 


( J _ 30 ’ 

\ J 20y 


sincu + 0 ljij (A21b) 


2a 2 (l - e 2 ) 2 


2 a 3 (l - e 2 Psin i 


(A 21 a) 


di _ c °s i 

“ 2a 3(i . ?2 f 


1-2 sin^l’lcOE (o + o(-22 
^ / \ J 20/ 


(A 21 c) 


de 

dt 

da 

dt 


3 A J 30 sin * 5 . 2 t\ , J J 30 

- — : =-(1 - t- sin 1 Jcos <n + 01 - — 

> V* 20/ 


2a 3 (l _ ||2) 2 ' * 


= 0 


( A21d) 
( A21e ) 


The change in the elements which are periodic in ft ’ can be obtained by 
considering the terms in equations (A9), which contain ft’. These equations 
can be put into an approximate equivalent form that can be integrated by 
replacing the value of i in the coefficients by its mean value taken over the 
period of the moon’s rotation about its polar axis. Neither e nor a need 
be considered as mean values since they do not have a change which is periodic 
in ft’ . 
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If the total rates of change in the elements (excluding short-period 
effects) as given by equations (A9) are separated into terms denoted by a sub- 
script 2, 3, or k according to whether they are periodic in 0,' , periodic in , 

or secular, they can be written in terms of the mean elements as follows: 


_ / da\ _ / da\ _ 
2 \dt /3 \dt/4 


= 0 


f 2s) = o 

idt /2 


(A22) 

(A23a) 


de 


5 nE m J 30 sin * 5 j 2 - i 

— — jll - ^ sin i )cos o> 


dt/3 2a5 _ - 2 )‘ 


(A23b) 


dt/4 


(A23c) 


-) - 

dt/2 


. ***** Sln 1 sin 2a' 
a 2 (l - e 2 ) 2 


(A24a) 


“I 

dt/3 


3 nJ 30 R m 


cos i 


2a 5 (l - e 2 ) 3 


1 - ^ sin 2 I^cos os 


(A24b ) 


dt/4 


(A24c) 



2 

JnJ 22^m / 

dt/2 

" a 2 (l - e 2 )^ ' 

W\ 


.dt/3 

2a3(l - 5 2 ) 3 

'doA 

3 aJ iO R rfl / 

,dt/4 

a2(l - e2) 2 ^ 


(l - | S in 2 i) 


cos 2Q 


A 72 2 t 


o> 


5 „j„2t 


(A25a) 


(A25b) 


(A25c) 
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p _ 

cos i 

a 2 (l - e^) 


cos 2£l x 



5nJ ? oR ” g - 2. sin^i)sin <o 

2aJ(l-e2) 3slnlV “ J 


(A26a) 


(A26b) 


/dn\ = _ 5PJ 20 Rg cos 1 

WA 2a 2(! _ S 2) 2 


The secular changes can be obtained from equations (A22), (A23c), (A24 c), 
(A25c), and (A26 c) as 


The rates periodic in 
sion 


(Aa-V = 

(Ad)i 4 . = 

(^1)4 = 

£2 T can 


(Ae)^ = (Ai)^ = 0 


3 nJ 2Q R m 


.2(1 _ 52) 2 


1 - r- sin 2 i Jt 

4 


3 nJ 20 R m cos 1 


(A27a) 

(A27b) 

(A27c) 


2a 2 (l - e 2 ) 2 

be integrated by use of the approximate expres- 


dfl* 

dt 


n l 


(A28) 


(where nj_ is the rate of rotation of the moon about its polar axis) to get 


(Ai) 2 = 

JnJpoR^ sin i , 

f 

cos 

= (cos 2£2 - 


2a 2 (l- e 2 ) nj_ 



(Ad) 2 = 

3c j 22 r^ / 

L - p- sin^iVsin 2£2 ! - sin 2Q { 

2a 2(l - e 2) 2 ni V 

4 

/ \ 

(Aft) 2 = 

2 — 

3nJ22 R m cos ^ 

/ . 

t \ 

2a2(l - e2)nT 

fsin 2£2 

- sin 2H 0 J 


(A29a) 


(A29b) 


(A29c) 
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The rates periodic in cd can he integrated by assuming that the major 
rate of change of <x> is due to its secular rate 


dn> 

at 



3nJ gQ R 2 / 
a 2 (l - e 2 ) 2 ^ 


sin 2 i 


*) 


This gives 


(Ae)j 


JjoR m sin 

2aJ20 


T 

1 


-(sin co 


sin cjo 0 ) 


(A30) 


(431a) 


J^oRm® cos i 

(Ai), = — ^ — 7 5 “t( sin co - sin co 0 ) (A31b) 

2 aJ 2 0 v .1 - e ) 


(Ao) 3 

(Aq)3 


j 3C)R] 


L m 


sin 2 i 


-2 2 ? 
e cos x 


2aJ 2 o(l " e 2 ) e sin I 


(cos CO - COS C0 o ) 


7 — Y cot i(cos co - cos co 0 ) 

2aJ 20 (l - e 2 ) 


(A31c) 


(A3ia) 


To obtain a comparison of the relative magnitude of these changes, the maximum 
amplitudes of the rates for a given orbit are given in the following chart: 


^ Subscript 

Element 

2 

3 

4 

de/dt, per day .... 
di/dt, deg/day .... 
cko/dt, deg/day .... 
cLft/dt, deg/day .... 

0 

0.048 

0.090 

0.124 

0.00125 

0.195 

0.278 

0.103 

0 

0 

1.142 

0.635 


For computation, the following orbital elements are considered: 

a = 2224 kilometers 
e = 0.1972 

I = 21° 
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J 20 = 2 -073 
J 22 = -2.03 ; 
J 30 = -9.3 X 


X 10-^ 

X 10-5 
10" 5 
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SECOND-ORDER ANALYSIS 


The following analysis is second order in the sense that an approximate 
expression which includes both secular and long-period (periodic in co) effects 
are used in the expression for cko/dt. This expression for dco/dt is then 
used in the determination of Arp rather than the approximate expression 

( — ] . According to equations (A25b) and (A25c), the sum of the secular and 

\ dtyi + 

second long-period change in the argument of pericentron is 


&JD 


3nJ 2Q B^ 


1 - y- sin 2 i 

k 


1 + 


J 30 Rm sin co sin 2T _ § 2, 


cos^I 


2J 2Q a(l - e 2 ) 


dt a 2 (l-e 2 ) 2 ' 

For convenience, the following definitions are presented: 


e sm l 


0)c = 


5 (l - 52) : 


S Pt 

1-7- sin i 


(bi) 


(B2) 


f = 


JjoRm 


2 - -2 2 - 
sin i - e cos i 


2J 20 a(l - e 2 ) e sin i 


Substituting equations (B2) and (Bj) into equation (Bl) gives 

dw • r ^ n \ 

— = cjd s (1 + f sxn a)) 
dt 


Integration of equation (b 4) yields 

r 


do) 


“S -CDo 

and relates t and <x> by the relation 


+ f sin cd 


(B5) 


(B4) 


(B5) 


t = 


CD, 


N5T 


f + tan §\ It + tan ^ 

rtB + tan - -*-! - tan“- M 


^7 


(b6) 
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where 8=1 when o> 0 < jt < CD and 8 = 0 for all other values of od and 
odo- In all cases cDq ^ cd < cdq + 2jt. 

To the second order in f, equation (B6) can be expanded to obtain 


t 



+ f(cos CD - 


cos CDq) 


+ 




(B7) 


Equation (B7) is accurate only when f is small compared with unity. Since 
the long-period change in the pericentron altitude is given as 


Arp = -a Ae 


(B8) 


an expression is needed for the long-period change in eccentricity. This 
expression is obtained from the rate de/dt given in equation (A21d) which 
can be written as 


de 

at 


sin 1 . 

CD- 

2aJ 20 S 


cos CD 


(B9) 


Substitution of dm/dt from equation (B^) into equation (B9) gives 

sin 


Ae = - 


2aJ20 

J^pgm sin i 


1 PCD 


COS CO 

1 + f sin cjl> 


dn) 


2J 2Qaf 


log 


1 + f sin cd 
1 + f sin ov 


(BIO) 


which gives the following long-period change in the pericentron altitude 



(Bll) 


It may be seen that equation (Bll) reduces to the first-order result for Arp 

whenever f is small. This condition is illustrated by the expansion of equa- 
tion (Bll) in a Taylor T s series, which is 
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_ J g0 R m sini r _ i f2sin ^ + i f3 sin 3 a) 

p 2J 20 f L 2 5 


- f sin a> 0 + i f 2 sin 2 cu 0 - j f^sin^a) 0 + C^f^) 


J^qRh! sin i 

Ar p = (sin to - sin cdq) 

2J20 


1 - —(sin co + sin cdq) 


— ^sin 2 oD + sin a> sin <n 0 + sin 2 ^ + o(f5) 


Caution should be used in the application of equation (B12) since it is 
restricted to small values of f. 
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TIME HISTORY OF PERICENTRON LOCATION AND SOLAR POSITION 


The dependence of latitude and longitude of pericentron, solar positions, 
and values of argument of pericentron on inclination can be obtained from the 
following sketch: 



Application of spherical trigonometry to the sketch results in the following 
relations : 


sin 0 
cos 0 
sin 0 


= sin i 
cos 0 = 
cos 0 = 


sin 0) 

-sin ft' sin a) cos i + cos ft* cos oo 
cos A* sin cd cos i + sin ft' cos oo f 


cos(cr 


0)cos 



J 


(Cl) 


Given initial values of co and fl ' , the initial latitude and longitude 
of the argument of pericentron as a function of i can be found from equa- 
tion (Cl) as 


sin 0o = sin i sin cd 0 

sin Aq sin cd 0 cos i + cos SIq cos o> 0 
cos 0 O = - — — 

^1 - sin 2 i sin^jDo > 

cos n' sin <n 0 cos i + sin cos co 0 
sin 0 O = — ■ - — 

1 - sin 2 i sin 2 o)o ^ 


(C2) 


The initial values of the solar positions for the initial, longitude and 
latitude can be obtained from the expression 
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cos(cr 0 - 0 O ) = 


2 v jl - sin 2 i sin 2 ^ 


(C3) 


The subsequent change in <f) with time is obtained by substituting the fol- 
lowing first-order expressions for cd into equation (Cl): 


CD = COq + 


3 ^ 20 ^ 

a^(l - e^)' 

The change of 0 with time is obtained as 




sin 2 i)t 


(c 4 ) 


e = e 0 - (13.2 - n)t = e 0 - 


13.2 + 


3nJ2QR; 


m 


cos 1 


2a‘ 


! (l - 


(05) 


The change in solar position a, with time is obtained by substituting into the 
last equation in equation (Cl) the time-dependent values of f and 9 just 
described. 
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